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C^*) ' We apply the technique of far-field interferometry to measure the properties of surface waves 

generated by two-dimensional (2D) single subwavelength slit-groove structures on gold films. The 
effective surface index of refraction n sur f measured for the surface wave propagating over a distance of 
■ more than 12 fim is determined to be n sur f = 1.016±0.004, to within experimental uncertainty close 

to the expected bound surface plasmon-polariton (SPP) value for a Au/Air interface of n spP = 1.018. 
O |i We compare these measurements to finite-difference-time-domain (FDTD) numerical simulations of 

the optical field transmission through these devices. We find excellent agreement between the 
measurements and the simulations for n sur f. The measurements also show that the surface wave 
, propagation parameter fc sur f exhibits transient behavior close to the slit, evolving smoothly from 

greater values asymptotically toward fc spp over the first 2-3 fim of slit-groove distance x ss . This 
behavior is confirmed by the FDTD simulations. 
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O ' !• INTRODUCTION 

o ' 

Recent measurements [1, 2] have characterized surface waves arising from optical excitation of a series of subwave- 
\^ ■ length slit-groove structures fabricated on silver films. The amplitude, wavelength and phase of these surface waves 
have been measured over a slit-groove distance of a few microns. After an initial rapid amplitude decrease extending 
. to ~ 3 /im from the slit edge, the interference fringe persists over several microns with a near-constant amplitude 
and a contrast ~ 0.3. The "near-zone" of rapidly changing surface wave character indicates initial transient phenom- 
ena, while the longer-range "far-zone" settling to constant amplitude and contrast is the signature of a guided mode 
surface plasmon polariton (SPP). The measured interference fringe wavelength A sur f = 814 ± 8 nm resulted in the 
| determination of an effective index of refraction n sur f = 1.047, significantly higher than expected from conventional 
Qh theory [3], n spp = 1.015. The question naturally arose as to whether this result is simply a consequence of interference 
fringe sampling over an interval predominantly in a transient near-zone peculiar to silver films, or was related to the 
specific surface properties of the silver films deposited on fused silica [4] , or may indicate a more generic phenomenon 
related to the transient properties of surface waves, generated by subwavelength slits, within 2-3 wavelengths of the 
edge of origin. Subsequent investigation of the physical-chemical properties of the silver films confirmed that they 
were free of surface contaminants, and FDTD simulations showed excellent agreement with the experimental results 
in silver films [5, 6]. 

In order to explore these surface waves further we have carried out a series of experiments similar to those already 
reported but using evaporated gold films instead of silver. We performed only "output-side" experiments [1] because 
they yield a phase modulated interference fringe less susceptible to noise than the amplitude modulation signal of 
the "input-side" experiments [2] . In addition we have compared the measurements to detailed field amplitude and 
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FIG. 1: Diagram showing interfering wavefronts, optical path difference between Et and E g , and far-field detection. Slit 
dimensions are 100 nm width, 20 /im length. Groove dimensions are 100 nm width, 100 nm nominal depth, 20 fim length. The 
evaporated gold layer deposited on a 1 mm fused silica substrate has a 400 nm nominal thickness. 
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FIG. 2: Goniometer setup for measuring far-field light intensity and angular distributions. See text for description. 

phase maps generated by FDTD numerical solutions [7, 8] to the Maxwell equations in the vicinity of the slit-groove 
structure. Both experiments and FDTD simulation show that the surface wave exhibits transient properties in 
wavelength and amplitude in the near-zone. This transient behavior can be interpreted in terms of surface modes all 
of which dissipate beyond the near-zone except for fc spp , the bound surface mode. 

II. EXPERIMENTAL SETUP 

The structures were fabricated with a focused ion beam (FIB) station as described previously in [2] . The experiments 
were performed on the same home-built goniometer setup used in the silver studies. For convenience, Figs. 1,2 of [1] 
are reproduced here to show the principle of the measurement and the schematic arrangement of the apparatus. 
Output from a temperature stabilized diode laser source, is modulated at 850 Hz by a mechanical chopper, injected 
into a monomode optical fibre, focused and linearly polarized (TM polarization, H-field parallel to the slit long axis) 
before impinging perpendicularly on the matrix of structures mounted in a sample holder. The beam waist diameter 
and confocal parameter of the illuminating source are 300 /zm and 33 cm, respectively. The sample holder itself is 
fixed to a precision x-y translator, and individual slit-groove structures of the 2-D matrix are successively positioned 
at the laser beam waist. A photodiode detector is mounted at the end of a 200 mm rigid arm that rotates about 
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FIG. 3: Left panel: Points are the measured fringe phase ip(x sg ) as a function of slit-groove distance a; sg . The straight-line fit 
is ipo = fc S urf Xs S + <fiint with constant slope k SUT f and intercept ipi nt . Gap in the data in left and right panels of this figure and 
in the left and right panels of Fig. 4 is due to defective structures in this interval. Right panel: Fringe amplitude rjo — 1/2C as 
a function of a; sg where C is the interference fringe contrast. 



the center of the sample holder. A stepper motor drives the arm at calibrated angular increments of 1.95 mrad per 
step, and the overall angular resolution of the goniometer is ~ 4 mrad. The photodetector output current passes to a 
lock-in amplifier referenced to the optical chopper wheel. Data are collected on a personal computer that also controls 
the goniometer drive. 



III. RESULTS AND ANALYSIS 



A. Measurements 



With the detector rotated perpendicular to the structure plane (6 = in Figs. 1, 2) the expression for the normalized 
detected intensity I / Io as a function of slit-groove distance x sg is given by 

— oc 1 + T)l + 2r] cos(k x x sg + ip int ) (1) 



Io 

where rjo is related to the fringe contrast C through 



c = ^r~2 ( 2 ) 

1 + To 



In the argument of the cosine term k x ~ 27r/A c ff relates the propagation parameter of the surface wave to the effective 
surface wavelength, and the "intrinsic phase" (p- mt is any phase contribution not directly due to the propagation path 
length x sg . It may be associated with phase shifts at the slit or groove structures. The interferometry measurements 
were carried out on 4 separate substrates, each substrate containing about 50 structures in which the the slit-groove 
separation was systematically varied from 50 nm to more than 12 /jm in increments of 50 nm. The left panel of Fig. 3 
plots the measured interference fringe phase against the slit-groove distance x sg . The fitted value for fc SU rf = 27r/A SU rf 
determines the effective surface index of refraction n sur f , and extrapolation to zero slit-groove distance a; sg determines 
the intrinsic phase (fii n t- The fit from the left panel of Fig. 3 yields 

n SUTi = = 1.016 ± 0.004 and (p int = 0.35tt ± O.OItt (3) 

^surf 

The right panel of Fig. 3 plots 770, the amplitude factor of the interference term in Eq. 1, as a function of slit-groove 
distance. From Eq. 2 this amplitude factor can be expressed in terms of the fringe visibility or contrast C as 

Vo= 1 ~ V ^ ^\c, C«l (4) 
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FIG. 4: Left panel: Fringe phase difference tp(x ss ) — ipo as a function of slit groove distance x sg . Deviation in the near-zone 
from Lpo indicates that early, transient fringe oscillation is slightly greater and approaches ipo asymptotically in the far-zone 
beyond ~ 2 fim slit-groove distance. Right panel: Same data as shown in left panel but on an expanded scale of slit-groove 
distances to emphasize the curvature in ip(x 3g ) — ipo in the near-zone. 



Although the fringe contrast shown in the right panel of Fig. 3 is about a factor of 5 below that measured for silver 
structures [1], the same rapid fall-off in the near-zone, x sg ~ — 3 /im, followed by a near-constant contrast beyond is 
observed. This contrast behavior is evidence of surface-wave transient phenomena in the near-zone. More evidence of 
this transient behavior is shown in left and right panels of Fig. 4 that plot tp(x) — (po vs. x sg , where ipo = k snT fX sg + ipi a t 
is the best- fit linear trace in the left panel of Fig. 3 over the range of slit-groove distances out to 12 /im. A pronounced 
departure from the asymptotic value of ipo is evident in the near-zone of slit groove distances, indicating that the 
fringe oscillation frequency is initially somewhat greater than the SPP value and smoothly decreases to it beyond the 
near-zone. 



B. Numerical simulations 



The time-dependent Maxwell equations are solved numerically using an FDTD non-conformal grid refinement 
method in Cartesian space coordinates. The methodology is described in greater detail in Refs [7-9]. 

Figure 5 shows a field map of the z-component of the electric field amplitude, and Fig. 6 shows a field map of the 
y-component of the magnetic field amplitude. The |E Z | map clearly shows the dipolar charge distribution concentrated 
at the corners of the slit on the input and output planes of the structure. These corner charge concentrations result 
from currents induced on the input side of the gold film by the magnetic field components |H y | shown in Fig. 6. The 
incident light propagates from below through the fused silica substrate onto the gold film and through the 100 nm 
wide, 400 nm thick slit. The incident light is TM polarized and the guided mode propagating along the ±z direction 
within the slit sets up a standing wave resulting in a high |E Z |, |H y | amplitudes at the output plane. The groove is 
at the output side of the gold film; and, in the simulations depicted, the distance between the center of the slit and 
the center of the groove is 3.18 /im. The absolute value of the z-component of the electric field amplitude |E Z | just 
above (or just below) the gold film is proportional to the surface charge density at the film surface. Note that on the 
output side surface, in the region between the groove and the slit, the surface wave excited at the left edge of the 
slit travels to the groove, is reflected from the groove's right edge, then interferes with itself. The standing wave is 
clearly visible on and near the output side surface. Within the slit, on the vertical walls, E z is fairly strong as well. 
Here, however, E z is parallel to the metallic surface, and its presence within the skin depth of the slit walls does not 
signify the existence of surface charges; instead, the E z -field in this region is responsible for the surface currents that 
carry the charges back and forth between the entrance and exit facets of the slit. Figure 6 shows the magnitude of 
the magnetic field |H y | over the same region as Fig. 5. Interference fringes between the incident and reflected beams 
on fused silica substrate are clearly visible. Note also the interference fringes between the excited evanescent waves 
(mainly SPP) and the incident beam on the entrance facet of the gold film adjacent the substrate. Inside the slit, 
|H y | shows a dark band; this is caused by interference between the upward-moving guided mode within the slit and 
the reflected, downward-travelling mode. 
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FIG. 5: FDTD simulations for slit-groove center-to-center distance of 3.18 /im, slit and groove widths 100 nm, groove depth 
100 nm and gold film thickness 400 nm. Map shows |E Z |, z-components (perpendicular to input and output facets) of the 
electric field amplitude in the vicinity of the input and output surfaces. 




x (|iim) 

FIG. 6: FDTD simulations for slit-groove center-to-center distance of 3.18 /im, slit and groove widths 100 nm, groove depth 
100 nm and gold film thickness 400 nm. Map shows |H y |, y-components (parallel to the slit and groove long axis) of the 
magnetic field amplitude in the vicinity of the input and output surfaces. 

In addition to the electric and magnetic field components E z ,H y at the surface, the light transmission efficiency 
in the z direction through the slit as a function of slit-groove distance x sg was calculated and is shown in Fig. 7. 
The transmission efficiency T is defined as the ratio of the z-component of the Poynting vector S° ut on the output 
side, integrated over x, to the total energy flux incident on the slit S™ . The red trace plots the transmission 

efficiency T = S° ut / S™ as a function of x sg at the output plane. As expected the transmission efficiency exhibits 
pronounced oscillations with a rapid decrease in amplitude in the near-zone followed by a constant amplitude oscillation 
out to 6 /im, the calculation limit. In the far-zone the blue trace fits the oscillations to a single cosine function, 
cos [2 (fcg^*f • x sg + <^j d )] . These oscillations, as can be seen in Figs. 5 and 6, arise from the superposition of waves 
launched from the slit and back reflected at the groove. Because these waves are counterpropagating along x, rather 
than copropagating along z, the intensity fringe frequency is twice the fringe frequency of the far- field interferometry 
results expressed by Eq. 1. This standing wave at the output plane results in 

T = ST 1 /ST (x 2 [1 + cos 2 (kgg ■ x sg + ^ t td )] = 2 { 1 + cos 2 [<p (x sg )} } (5) 

Taking into account this factor of two in the argument of the cos term, the fringe oscillation in the asymptotic far- 
zone is in good agreement with measurement and the expected SPP. In the near-zone the oscillation tp (x sg ) exhibits 
a definite "chirp," and Fig. 8 plots the deviation from the asymptotic value tpo as function of x sg . Comparing Fig. 8 
to Fig. 4 we see that the FDTD results accord well with deviations in the interferometric fringes measured in the far- 
field. 
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FIG. 7: FDTD calculations of the transmission efficiency T — S° ut /Sl n as a function of a; sg . Red curve traces T, and the 
blue curve traces a cos (2fc^*f • x sg + ^inJ d ) fit to the oscillation in the asymptotic region. Note the decreasing transmission 
amplitude in the near-zone close to the slit edge and the higher oscillation frequency compared to the asymptotic harmonic 



wave. Best-fit values for A^'f — 27r/fc^*f = 839 nm and intrinsic phase </?fnt d = 0.557T rad. 
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FIG. 8: Phase difference <p(x sg ) — <po as a function of x as , analogous to the right panel of Fig. 4 but derived from the FDTD 
simulation data. Residual "high frequency" oscillations in the phase difference, believed to be due to numeric artifacts in the 
FDTD results, have been smoothed. 



IV. DISCUSSION 



Surface wave in the far-zone 



The index of refraction for the bound surface-plasmon-polariton n spp is given by the Raether formula [3] 



fc spp 



e m ed 



Ed 



(6) 



where e m and £d are the real parts of the dielectric constants of metal and dielectric at the interface on which the 
surface wave propagates. Interpolation of reflectivity data for gold [10] at 852 nm yields eau = —28.82 and from Eq. 6 
the surface index of refraction for the surface plasmon polariton at the gold-air interface is e spp = 1.018. The measured 
surface index of refraction reported here, to within experimental uncertainty and in the far-zone, is in accord with 
the SPP prediction. The results from the FDTD calculations are also in agreement with the experimental results 
and the SPP predictions. It appears therefore that in the far-zone, for both silver/air and gold/air surfaces, far-field 
interferometry and FDTD calculations show that the surviving long-range surface wave is indeed the expected bound 
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surface plasmon polariton. Table I summarizes the relevant parameters, far-field interferometric measurements, and 
finitc-diffcrence-time-dependcnt (FDTD) numerical simulations for gold and silver. 



TABLE I: Summary of A sur f, n SUI { and n spp determined from far-field interferometric studies and FDTD simulations in gold 
and silver 



A sur f Tlsutf 72spp 

Au,Ref. [11] -28.32 839 ± 6 nm 1.016 ±0.004 1.018 

Au,Ref. [12] -31.62 839 nm 1.016 1.016 

Ag,Ref. [13] -33.27 819 ± 8 nm 1.04 ±0.01 1.015 

Ag,Ref. [14] -33.27 814 ± 8 nm 1.05 ±0.01 1.015 

Ag,Ref. [15] -33.98 837 nm 1.017 1.015 



B. Surface wave in the near-zone 



In the near-zone both experiment and numerical simulation show that the surface wave deviates from pure SPP 
behavior. The effective propagation parameter k Xl originating near the slit edge, appears greater than k spp and 
evolves smoothly to the bound mode over the near-zone interval of a few microns. This behavior may be interpreted 
either in terms of initial excitation of a composite evanescent surface "wave packet" in fc-space at the slit edge [16], 
followed by subsequent decay of all surface modes except the bound fc spp mode or in terms of detailed field matching 
at the boundaries within the slit and near the slit edges [17]. These two points of view both invoke evanescent modes 
k x > k S pp m order to satisfy boundary conditions in the vicinity of the slit edge, but standard wave-guide theory 
dictates that only the SPP mode is stable against phonon coupling to the bulk metal or to radiative decay. We can 
estimate the the surface distance over which the dissipation occurs by appealing to the standard Drude model of a 
metal that expresses the frequency dependence of the dielectric constant e(lu) in terms of the bulk plasmon resonance 
uj p and a damping constant T. 

In Eq. 7 e is the permittivity of free space and = e(oj — > oo) is the dimensionless infinite frequency limit of the 
dielectric constant. The Drude model is based on a damped harmonic oscillator model of an electron plasma in which 
the electrons oscillate about positive ion centers with characteristic frequency w p , subject to a phenomcnological 
damping rate T, normally assumed to be due to electron-phonon coupling. Values for T are typically ~ 10 s , and 
in fact for Au the value is 1.02 x 10 14 s _1 . For a wave propagating on the surface with group velocity ~ 3 x 10 8 ms _1 
the expected decay length ~ 3 /jm, consistent with the measurements and FDTD simulations. 

In summary the picture that emerges from far-field intcrfcromctry and FDTD simulation studies of these simple 
slit-groove structures on silver and gold films is that in the near-zone of slit-groove distances, on the order of a few 
wavelengths, the surface wave consists of a composite of several evanescent modes all of which dissipate within this 
near-zone. Only the bound, stable SPP mode survives into the far-zone, and in the studies reported here we have 
observed essentially constant SPP amplitude out to 12 /im, a distance limited only by our fabricated structures. 
Earlier measurements [2] indicate that absorption and surface roughness scattering should permit propagation lengths 
as far as ~ 80 — 100 /im. 
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